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1> ' Abstract 
<N . 

. We obtain a family of first-order symmetric hyperbolic systems for the Bianchi 

I equations. They have only physical characteristics: the light cone and timelike 

' hypersurfaces. In the proof of the hyperbolicity, new positivity properties of the 

00 , Bel tensor are used. 

PACS number: 04.20.Ex 

u 

1 Introduction 

^ I First-order symmetric hyperbolic (FOSH) systems have recently been the subject of 

■ intensive research, not only for their theoretical interest in the Cauchy problem, but 

also for their applications to numerical relativity [0. Such systems constitute a pow- 
erful procedure that can be used to establish a well-posed initial-value formulation for 
physical systems, since most of them can be cast into this FOSH form and, within 
this framework, theorems of existence and uniqueness of solutions are applicable (see 
also [0). More particularly, the Bianchi identities (investigated for a long time mainly 
for its usefulness in the study of the gravitational field propagation) are the basis from 
where some FOSH systems have been extracted (see and references therein). Many 
of these results have been obtained by making use of the vacuum Bianchi identities, 
however the non- vacuum ones can also be applied to the construction of FOSH systems 
for the Einstein field equations 0. 

In this paper a family of FOSH systems for the Bianchi equations is constructed. 
They may be useful in deriving new hyperbolic systems for Einstein's field equations, 
analogously to what is done in 0], where Bianchi identities are linked to equations for 
the spatial metric and connection. We begin in section 2 by considering the Bianchi 
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equations applied to a double 2-form which does not satisfy a priori the additional 
Riemann symmetries. From these equations we extract a FOSH system of evolution 
equations having no unphysical characteristics. Actually, their only characteristics are 
the null cone and timelike hypersurfaces. In section 3 we briefly investigate the set of 
constraints that completes the hyperbolic system found in section 2. Both together, 
the evolution equations and the constraints, are equivalent to the original Bianchi 
equations. Finally, in section 4 we show that the full Riemann symmetries (together 
with the Einstein field equations) are conserved under the evolution. In the appendix 
we prove new positivity properties of the Bel tensor that are needed for the proof of 
the hyperbolicity. 



2 Hyperbolizations of Bianchi equations 

We recall that the Bianchi identities for the Riemann tensor read: 

V[y-Ra/3]A^ = , (1) 

from where it follows (by contracting, using the Riemann symmetries and the Einstein 
field equations with SnG = c = 1): 

^^RapX^i = 2V[^ (T\]a — -Tg\]o^ = JaXfi , (2) 

Ta/3 being the energy- momentum tensor. For the system (|l|,||), called the superior 
order field equations by Lichnerowicz, Einstein's field equations can be considered as 
initial data if they hold at the initial spatial hypersurface S, they also hold in a 
neighborhood of S. In the last section we will see that the symmetries can also be 
taken as initial data. 

Consider now the set of equations (||,0), which will be called Bianchi equations in 
what follows, applied to a double 2-form Kapx^j. = -^[«/3][a^(]- We do not suppose it to 
have the additional Riemann symmetries K^aiBXjfi = at the present moment. The 
Bianchi equations for Kapxn are then: 

'^^KajSXfi = JaX/i , (3) 
V[,/C/3]A^ = . (4) 

This is an overdetermined linear first-order system of 48 equations for 36 unknowns 
which is not manifestly symmetric hyperbolic. To handle more easily equations (^^, 
it is better to write them in the following compact form (in a similar way to what is 
done in p| for other systems): 

Ar''''V,K,s.p=ji, (5) 
where the index / lives in the space of equations, / = (a[A/i], [Ayu]), and 
^e,s., ^ {6^l'6l%5i^,6lJ^:syf,6^^) , jj = (J..„0) . 
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A tensor j's'^^'p' jg gg^i^j hyperbolization of if the matrix W^''^'p'}{i^^p} = 

^7y5V'p'^^75-P acting 

over the space of double 2-forms satisfies the following prop- 
erties: i) it is symmetric, and ii) there exists a 1-form such that n^Q^ is positive- 
definite. By means of such a hyperbolization, we can transform (§) into: 



•ySap 



Jl 



which is now a FOSH system and, therefore, it admits the theorems of existence and 
uniqueness of solutions. 

The set of 1-forms satisfying property ii) has a clear physical meaning (see 
for instance) related to the propagation of the physical fields involved and to the causal 
character of the system. In fact, a FOSH system is said to be causal if every satisfying 
ii) is a. future-directed 1-form. 

For the Bianchi equations (13) the following family of tensors are hyperbolizations: 



f'S'a'p' 



.a[yf]j^M[a'p'] _S_^[u.a\[Y.5']\f^]j^M[a'p' 

' 2 



y^Xpa'p' 
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uiPyP') _ 



where t, u and v are arbitrary timelike future-directed vectors. To check property i) 
we compute, for any two double 2-forms S^sap and T^sap, the quantity 



{7'<5'cr'p'}{7'5o-p} q , ,T , 



f qa 13 rp I rpa (3 q 



qa rp I rpci P q 



crp 



:9vp 



QCtPy rp I rpaPy Q 



^Oa/S-yp 



f q a/S-yrp I rp a(3y q \ I 

[Jap \ ^y pap-y ~r -t ^Jpaf^y j i 
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from where the symmetry is clear. Property ii) is demonstrated in the appendix by 
proving that n^Q*^ is a positive-definite quadratic form for any timelike future- 
directed 1-form. This also shows that the system is causal. 

Let us now study the characteristics of the hyperbolized Bianchi equations found 
above, which finally read: 



(6) 



Recall that a 1-form (a is, by definition, characteristic for the system (P) if det(CeQ'') = 
0. To study the kernel of the characteristic matrix, first notice that h[ai3][xp] (considered 
as a matrix acting over 2-forms) has non-zero determinant, that is: 



h 



al3][Xp] 



(7) 
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This can be easily proven by contracting (|^ successively with u and v. Therefore, 

The non-trivial solutions of this system are found to be of two types: i) CaC" = ^^id 
a) Cq^" = 0. In the first case, the solutions are the double 2- forms K^px^ satisfying 
C"i^Q,/3A^ = and C,[uKai3]\^ = 0. In the second one the solutions are those K^px^ 
having in its two first indices the structure A{t /\ Q) + B * {t f\ Q) ( "*" is the usual dual 
operator). Therefore, as remarked above, has only physical characteristics: i) the 
light cone and ii) timelike hypersurfaces. 

3 Equivalence with the original system 

From the Bianchi equations we have constructed the 36 propagation equations (H) (in 
the sense that they contain derivatives in any timelike direction n^). In this section 
we will show that these equations, when completed with the constraints, are equiva- 
lent to the original system Taking n^j as the normal to an initial- value spatial 
hypersurface S and r;^"^*^ the spacetime volume 4-form, the constraints are: 

n"(v^K,/3A^- J^Am) =0, (8) 
n,r;^"^'^V,i^„^A;. = 0, (9) 

where all the derivatives appearing in these expressions are tangent to S. This set 
of constraints is complete (there are 12 equations) and also integrable provided that 
'^"'JaXfi = 0, which holds because of the Einstein field equations. 
Taking into account (0), equation (P) immediately implies 

If we multiply the previous expression by ts, we obtain {y^K^^px^ — J^xii) = 0, 
where is the orthogonal projector to ta- Taking into account (|) and the fact that 
tan"" 7^ 0, we obtain equation (|^). Equation (^ is achieved in a similar way: we first 
multiply (|l3) by T]^sapt^ obtaining P^ri'^^^^V eKapxiJ. = 0, and then we follow the same 
reasoning as before, but now making use of constraint (|^). This shows that the set of 
equations (||J|J^) is equivalent to (0,§). 

4 Riemann symmetries are conserved 

Until now, we have only dealt with double 2- forms, but, actually, we are interested 
in the Riemann tensor. Our purpose now is to show that Bianchi equations preserve 
the full Riemann symmetries. The complete set of Riemann symmetries is KajjXfj. = 
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K[ai3][Xfi] and K[ai3x]^ = 0. In order to construct an appropriate system of differential 
equations, let us set the following definitions: 

where K^is = K^^^^. Then we extract, from equations (0,^) and the energy-momentum 
conservation equation V^T^p = 0, the following system for the previous defined vari- 
ables Ea/BXf, and Tj^/b- 

'^""^apXfj. + 2 V[/3Sa]/, = , 

V^T.^p - iv„S = . (11) 

These are 32 equations for 32 unknowns forming a linear first-order homogeneous sys- 
tem whose only characteristic is the null cone. Now, we suppose that initially the 
double 2-form Ka/3x^ satisfies the full Riemann symmetries (Sq^a/x = 0) and that, in 
addition, satisfies Einstein's equations in the sense S^^ = 0. These initial condi- 
tions have no other solution but H^^a^j = ^ai3 = (see 0], for instance) and, therefore, 
the Riemann symmetries, likewise the Einstein field equations, are preserved by the 
Bianchi equations. 
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A Positivity property of Bel tensor 

To prove that n^Q^ is a positive-definite quadratic form we will make use of the standard 
spinor calculus First, taking into account the symmetries of Kaisxfi, we decompose 
it as follows: 

Kupap = XNMSR^N'M'^S'R' + NMS'R'^N'M'^SR 
+ Xn'M'S'R'^NM^SR + ^ N'M'SR^NM^S'R' i 

where 

XNMSR = XiNM){SR) ) "^NMS'R' = '^{NM){S'R') ■ 
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Now we proceed to calculate 

-]^g,,K'^^\K^p,, - ^g^pRJ^P-'K^^p, + ^g^.g^pK^^^sK"^^') n^t^u'^vP , (12) 

which is the Bel tensor for the double 2-form K^p^s contracted with four timelike 
future-directed vectors, n, t, u and v. In fact, this is a little more general since we have 
not required symmetry between indices af5 and A/U. In order to express ([T2|) in spinor 
form it is more convenient to rewrite it in the following form: 

and now, the dual terms appearing in this expression have simple spinor expressions 
(see [^). Performing the calculations we finally obtain 

= 4 {xnmsrXn'M'S'R' + '^nms'R''^n'm'sr) n^fu'^v'' , (13) 

which is always positive, and zero if and only if Xnmsr = ^^nms'R' = 0. This is proven 
by using usual methods (see, for instance, p], |10[): let us consider a spinor basis (t^, oa) 
and its two associated null future-directed vectors, la = i-aIa' and = oaOa'- Given 
a timelike future-pointing vector z and a future-directed null vector x, there exists 
another future-directed null vector y such that: 

z = ax + y = axAXA' + VaVa' , a > 0. 

Applying this decomposition to n, t, u and -u, ([13|) becomes a sum of positive terms, so 
the positivity is shown. Suppose now that the expression (0) is zero. We can decom- 
pose the four previous vectors in many differents ways by choosing their corresponding 
vector a; to be / or k. Then, for every choice we do, among the positive terms of the 
sum we will get a component of Xnmsr and "i^NMS'R' in the spinor basis mentioned 
above. Therefore, making all the possible choices, we will obtain all the components 
zero, which means that Kapx^ = 0. 
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